Dislocation Glasses: Ageing during Relaxation and Coarsening 
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The dynamics of dislocations is reported to exhibit a range of glassy properties. We study numer- 
ically various versions of 2d edge dislocation systems, in the absence of externally applied stress. 
Two types of glassy behavior are identified (i) dislocations gliding along randomly placed, but fixed, 
axes exhibit relaxation to their spatially disordered stable state; (ii) if both climb and annihilation 
is allowed, irregular cellular structures can form on a growing length scale before all dislocations 
annihilate. In all cases both the correlation function and the diffusion coefficient are found to exhibit 
ageing. Relaxation in case (i) is a slow power law, furthermore, in the transient process (ii) the 
dynamical exponent z ~ 6, i.e., the cellular structure coarsens relatively slowly. 

PACS numbers: 64.70.Pf, 61.20.Lc, 81.05.Kf, 61.72.Bb 



Since dislocations arc ubiquitous in crystals and they 
can have a dramatic effect on physical properties of ma- 
terials, understanding the morphology and dynamics of 
dislocations is of fundamental importance. They are the 
carriers of plastic deformation in solids [l| and are unwel- 
come impurities in single-crystal growth [3|. In 2d they 
play a crucial role in melting [3|, appear in Wigner crys- 
tals [J], and dislocations in vortex lattices [5[ determine 
properties of superconducting films [6| . Further examples 
where dislocations arise include magnetic bubble struc- 
tures [7|, charge density waves Q, colloidal lattices 0, 
and dusty plasma [lO|. 

The morphology of dislocations is one of the most stud- 
ied and least understood aspects of crystal defects. While 
there is a multitude of dislocation patterns, one can dis- 
cern a few generic types. Such are fractals and ladders 
in the presence of external shear [11|, as well as cellu- 
lar and diffuse configurations in unsheared crystals. The 
latter ones were discussed in the remarkable recent ex- 
perimental papers by Rudolph et al. [2, II4I with the con- 
clusion that climb mobility of dislocations has an impor- 
tant effect on morphology. Specifically, when there is 
climb then cells of some characteristic length scale can 
form, or, if climb is suppressed, dislocations freeze into 
a diffuse-looking random distribution. A recent work by 
the present authors corroborated this finding on the basis 
of mesoscale simulations [iSl . 

Just a glance at disordered dislocation configurations 
from experiments [3| and simulations [13|, ll4| suffices to 
raise the question, to what extent dislocations exhibit 
glassy properties. Indications of slow equilibration came 
recently from numerical studies of relaxing 2d disloca- 
tion systems [15|. A physical argument for the glass 
analogy is that the long-range interaction between dis- 
locations can change sign as function of relative angles, 
thus making frustration possible. Frustrated interactions 
are known to give rise to spin glasses in the presence of 
quenched disorder [16[. Furthermore, there is an exter- 



nal element of randomness in dislocation systems that 
arises from initial conditions. This reinforces the anal- 
ogy, because random initial conditions due to quench- 
ing from the high-T state arc an essential ingredient in 
structural-glass-formers, like in Lennard- Jones systems 
[l7|, and lead to slow growth of irregular domains even 
in systems as simple as the Ising model [18| . 

There is a growing body of knowledge about the dy- 
namics of glasses (see [18[ and rcfs. therein). A central 
feature of glassy dynamics is the dramatic slowing down 
while the temperature is lowered [19[ . Not only does the 
dynamics slow down, but it also exhibits ageing, widely 
considered as a hallmark of glassy dynamics [20| . There 
are several qualitatively different types of systems ex- 
hibiting ageing, such as (i) spin glasses with quenched 
disordered interactions, (ii) systems approaching a non- 
glassy equilibrium state, while random initial condition, 
possibly augmented with kinetic constraints, slow down 
relaxation, and (iii) structural glasses wherein long-time 
disorder originates in initial condition. 

In this Letter we study the dynamics of edge disloca- 
tions in 2d by mesoscale simulations [l3| , wherein disloca- 
tions are the particles interacting by their stress fields and 
obeying overdamped dynamics. While this is a strong 
simplification of real 3d dislocation networks, for FCC 
crystals under certain conditions it is a reasonable ap- 
proximation, and it is of course relevant for defects in 
2d lattices. Two main situations are considered, in the 
first one dislocations only glide along the axis of their 
Burgers vector, in the second one they can also climb 
perpendicular to it, and here we also allow annihilation. 
Climb has a dramatic effect, without it dislocations relax 
to a (meta)stable state, while climb with annihilation re- 
sults in an increasingly dilute system and eventually all 
dislocations vanish, corresponding to the dipolized state 
predicted long time ago [3|. In both situations we study 
the cases of a single and three glide axes. The effect 
of temperature on ageing is tested by the addition of a 



Langevin force. 

Our main observations are as foUows. (1) Disloca- 
tion configurations after a fast initial transient depend 
on whether climb is possible and on the number of glide 
axes. With climb, walls form for both the 1- and 3-slip 
systems, in the latter leading to a cellular structure with 
a pronounced length scale, whereas without climb typical 
configurations appear more disordered. (2) All systems 
exhibit ageing, that is, both the correlation function and 
the effective diffusion cocfScicnt of dislocations depend 
on the waiting time. (3) In the absence of climb the cor- 
relation functions decay by a power law while a stable 
state is approached. (4) Sufficiently high temperature 
restores non-ageing dynamics. (5) With climb, the time 
exponent z is identified, characterizing the growth of the 
mean domain diameter. 

First we discuss the details of the simulation. Three 
different possible Burgers vectors defined by the direc- 
tions ±(cos(TO7r/3),sin(TO7r/3)), to = 0,1,2 in a square 
area of size L x L was considered with periodic boundary 
conditions. This geometry emulates a hexagonal under- 
lying lattice. In order to calculate the stress field of a 
dislocation with Burgers vector (6a;, 5j,), we use the for- 
malism developed by Kroner [l|]. The stress function x(r) 
for this problem fulfills the biharmonic equation 



components of the ith dislocation are 



v^x 



C(h.,dy - hyd.,)8{x)6{y). 



(1) 



with C = 2/i/(l - v) in 3d and C = 2/i(l + v) in 2d, 
where [i is the shear modulus and v Poisson's ratio. From 
x(r) the stress components are a^x ~ ~'^v'^vXi '^yy — 
—dxdxXi o'xj/ = dxdyX- The force F per unit length 
acting on a dislocation with Burgers vector b is given 
by the Peach-Koehler formula F ~ {a b) x I in which 
I = (0, 0, 1) and the stress a is due to all other dislo- 
cations. Equation ^, discretized to M^ points, with 
periodic boundary conditions, can be solved for x by us- 
ing the Fourier transform (FT) of the partial difference 
operator, D„ = ^ sin (^) , where m = 0, . . . , M - 1. 
Then one obtains the FT of the stress function as 



Xi„^ = C{bxD„, - byDi)/[Df + Dl^W 
yielding the FT of the stress components dxx 



(2) 
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DiD„,xi^. In- 
verse FT gives the stress components at the grid points 
for a single dislocation. For the effect of all dislocations, 
in xim we replace b by the FT of the discretized Burg- 
ers vector field. While we keep track of the dislocation 
positions to high precision, the stress field is discretized 
to M = 256. Note that this recipe suppresses interac- 
tion between dislocations in the same box, but in our 
case such concurrence is rare. By applying fast FT, for 
up to 10000 dislocations the above algorithm is about 10 
times faster than the direct calculation of the pair inter- 
actions, for more details see [13|. Assuming overdamped 
dislocation dynamics, the v'i glide and v* climb velocity 



g,c 



= r,,,(F'a^_ J + y27TVC, 



S.c 



(3) 



where F^ is the Peach-Koehler force given above, Tg and 
Tc are the glide and the climb mobilities, respectively, 
o,g = bi/\b\i, and a^ is a unit vector perpendicular to a^. 
Thermal fluctuations at temperature T are included as 
white noise ^g^c with {^g,c{t)^g,c{t')) = S{t — t'). (So as 
to eliminate material parameters the dimensionless vari- 
ables r =^ r/L, t ^ tTgCh'^/ATiL'^, and T ^ AnT/Cb'^ are 
used henceforth.) 

In order to study the morphology of dislocations, we 
considered dynamics without and with climb, and in ei- 
ther case we took 1 as well as 3 glide axes. Simulation 
started out from random initial conditions of 10000 dislo- 
cations, specifically, axes were placed independently, with 
a uniform distribution, and the sign of the Burgers vector 
was also randomly chosen such that for each axis direc- 
tion they compensated each other. No shear is imposed 
through external traction, as we intend to study dislo- 
cations originating from some internal inhomogeneities. 
If climb is present, dislocations can get arbitrarily close 
to each other, so we annihilated dislocations of opposite 
Burgers vectors in the same box. Figures [T^-d show typ- 
ical configurations, without climb (a,b) and with climb 
and annihilation (c,d), in each case for 1 and 3 slip axes, 
respectively. Configurations (a,b) are near equilibrium, 
they show diffuse distributions, but in the single slip case 
the formation of walls is apparent, in accordance with re- 
cent analytic results 2l|. On (C;d) transient states are 
snapshot, where the predominance of walls is apparent. 
Dislocations eventually annihilate, in agreement with the 
classic prediction [3|] of dipolized equilibrium well below 
the melting point. As we have reported in I13|, walls 
form cells in the 3-slip case as seen in Fig.[Tli. The picture 
shows a close resemblance to some dislocation patterns in 
2d dusty plasmas [10| , highlighting the experimental rel- 
evance of our simulation. Furthermore, our results are in 
qualitative concordance with recent experimental studies 
[12| . where, in off stoichiometric GaAs, cell formation was 
attributed to increased climb mobility. So we can con- 
clude that first and foremost the presence of climb, but 
also the number of slip axes, strongly affect the morphol- 
ogy of dislocations. It is common for all cases presented 
here, that random initial conditions result in disordered 
configurations, calling for studies into glassy properties 
of dislocation systems. 

Aiming at quantifying the glassy character of disloca- 
tion dynamics we focus on ageing (see [18|). In particu- 
lar, we study the dependence on the waiting time of the 
correlation, or overlap, function 

C(iw, t) = ]^ I^exp [-\r,{t + tw) - r,;(iw)|/ro] , (4) 




FIG. 1: Snapshots of relaxed dislocation configurations with- 
out climb (Fc = 0): (a) single slip (b) 3-slip, and typical wall 
structures with climb (Fc/Ff, = 0.1): (c) single slip, (d) 3-slip. 



where somewhat arbitrarily ro — 5/^ 
tivc diffusion coefficient 
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|r»(t + tw)-r,(tw)|^ 



(5) 



where N^ is the number of dislocations remaining at the 
final observation time. Here C and D should be conceived 
as time t dependent quantities, recorded after a waiting 
time iw. 

First we consider the situation without climb, with 3 
slip axes, at zero temperature, and plot the correlation 
function on Fig. [2l The asymptotic value Coc{tw) = 
C(t„, oo) markedly depends on t^ (left inset). Generi- 
cally, beyond the inflexion point correlations decay ex- 
ponentially, followed by a power form t~^. Thus it is 
natural to consider the curves 1 — (1 — C)/(l — Coo) de- 
caying in time t from 1 to 0, which are found to collapse 
approximately to a master curve as function of t/tj . 
This is shown on the right inset, where the straight line 
indicates power law with (3 w 0.54. Note that the fact 
that tw and t have different exponents contradicts stan- 
dard scaling formulas in ageing (see [l8| ) , an anomaly for 
which we do not have an explanation and requires fur- 
ther studies. It is worth recalling here th at p ower decay 
of the elastic energy was found before in [l5| . So we can 
conclude here that there is a strong analogy with spin 
glasses inasmuch that there is a quenched randomness in 
the placement of the slip axes, which confine dislocations 
along the entire dynamics, and the relaxation towards a 
disordered ground state exhibits ageing. Next we study 
the dynamics of the 3 slip system without climb in the 
presence of a small Langevin noise, T = 0.025. Then the 




FIG. 2: Correlation function C{t^,t) for 3-slip, without climb 
at r = as function of t for different tw = 1, 2, ..., 2048. 
Increasing t^ leads to slower decay. Left and right inset show 
the asymptote and the scaling form, resp. 



system never fully relaxes, but before some stationary 
state is reached it also exhibits ageing, as demonstrated 
in Fig. [31 Plot (a) shows the correlation function de- 
caying cxpectedly to zero (we did not reach the asymp- 
tote due to limitations in computational power) for a se- 
quence of waiting times. On (b) the diffusion coefficient is 
shown, where lower lying curves correspond to increasing 
tw. Asymptotically we found D{tw,t) ~ Dq + Di{tv,)t~'' 
with Do ~ 2.210~^,7 w 0.8, see upper right log-log in- 
set. The power decay to Dq can be considered as a man- 
ifestation of slowed-down dynamics, which may appear 
for intermediate t-s as effectively subdiffusive. To test 
whether the remanent diffusion constant Dq is physically 
relevant or only due to finite size effects, we halved the 
effective linear system size by running a simulation with 
A'd = 2500. The lower left inset shows that the asymp- 
totic Dq docs not change significantly for varying system 
size, so ordinary diffusion indeed prevails. 

Finally, we considered dislocations with climb and an- 
nihilation for 3 slip axes. Here the number of dislocations 
A^d is a decreasing function of time, and after an initial 
transient the dislocations order into cellular structures, 
whose characteristic length scale grows while dislocations 
annihilate. Figure [3] shows the ageing effect in the over- 
lap function and the diffusion coefficient, different curves 
corresponding to different values of the waiting time. A 
small thermal force was included, presumably having no 
significant effect on dislocations in cell walls. The inset 
shows the correlation function at a sufficiently high tem- 
perature that suppressed ageing. As we neglected the cre- 
ation of dislocations, the process we followed is essentially 
a transient before all dislocations annihilate. Further- 
more, due to climb the dynamics is not constrained to the 
randomly placed glide axes, so disorder in the morphol- 
ogy is essentially due to random initial conditions. This 
situation is reminiscent to the transient observed in sys- 
tems with non-disordered equilibrium states but showing 
a transient with irregular domain coarsening if quenched 
from a sufficiently high temperature. While systems with 
interrupted ageing arc usually not considered as glasses. 
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FIG. 3: Correlation function C{t^,t) (a) and diffusion co- 
efficient D{tvi,t) (b) for 3-slip, witfiout climb, T = 0.025, as 
function of t for various iw-s. Upper inset in (b) shows the 
decay to the asymptotic Do = 2.2 fO~^ and lower inset gives 
D(0, t) for two system sizes. 
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FIG. 4: Correlation function (a) and diffusion coefficient (b) 
for 3-slip, with climb (Fc/Fg = 0.1) and annihilation, other- 
wise like in Fig. |3l Inset in (a) is the plot of the correlation 
function at a sufficiently high temperature showing that the 
dependence on the waiting time tw is practically lost. 



their transients do show ageing and, on short time scales, 
further glassy characteristics were detected recently [22] ■ 
A common property of those systems is the asymptotic 
scaling of the characteristic domain size i with time as 
I oc t^/^. For instance, z = 2,3 corresponds to the 3d 
Ising model with Glauber and Kawasaki dynamics, and 
the logarithmic growth in the random field model can 
be associated with z = oo, see [l^. Earher we found 
numerically that the number of dislocations decreased 
app roximately like A^d oc t^ with 5 approximately 1/3 
[13| . Hence, based on the scaling of all lengths by the 
dislocation number we can conclude that the dynamical 
exponent here is in the order of z « 6(1), where the er- 
ror estimate is subjective. This relatively large exponent 
means that the glassy transient lasts for comparatively 



long times and may explain why the cellular structure 
can persist under experimental conditions. 

We can summarize the physical picture of disordered 
dislocation systems as follows. In the absence of climb, 
where the external disorder is in the random but fixed 
placement of glide axes, there is an analogy with spin 
glasses in that the dislocation system has a multitude of 
near ground states, to one of which it relaxes. No in- 
creasing length scale was discerned at this stage of simu- 
lation. On the other hand, with climb, the ever growing 
cellular structure resembles the domain growth in, say, 
the Ising model launched after a quench from a param- 
agnetic state and showing glass-like dynamics. In each 
case there is a pronounced ageing effect in both the cor- 
relation function and the diffusion coefficient at zero and 
small temperatures. Interestingly, we did not encounter 
dislocation states analogous to structural glasses, proba- 
bly because dislocations moving without constraints can 
annihilate. Our results are, however, strongly suggestive 
that we can speak about dislocation glasses, where the 
types of glasses are distinguished by different morpholo- 
gies. 
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